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SWEEPING PROCESS BY PROX-REGULAR SETS IN 
RIEMANNIAN HILBERT MANIFOLDS 

by 

Frederic Bernicot Sz Juliette Venel 


Abstract. — In this paper, we deal with sweeping processes on (possibly infinite¬ 
dimensional) Riemannian Hilbert manifolds. We extend the useful notions (proximal 
normal cone, prox-regularity) already defined in the setting of a Hilbert space to the 
framework of such manifolds. Especially we introduce the concept of local prox-regularity 
of a closed subset in accordance with the geometrical features of the ambient manifold 
and we check that this regularity implies a property of hypomonotonicity for the proximal 
normal cone. Moreover we show that the metric projection onto a locally prox-regular 
set is single-valued in its neighborhood. Then under some assumptions, we prove the 
well-posedness of perturbed sweeping processes by locally prox-regular sets. 


1. Introduction 

1.1. A brief review of results on sweeping processes. — The aim of this paper 
is to extend the study of so-called sweeping process in the setting of a (possibly infinite 
dimensional) Riemannian manifold. 

Sweeping processes are specific differential inclusions of first order involving proximal 
normal cones to a moving set. They were extensively studied for the last years in the 
Euclidean space and more generally in a Hilbert space. 

More precisely, consider a Hilbert space H and a moving set t C'(t) on a time-interval 
I := [0,r] and assume that the set-valued map C : I ^ H takes nonempty closed values. 
A function u : / —)■ R is a solution of a problem of perturbed sweeping process if it satisfies 
the following differential inclusion: 

-h N(C'(t), u(t)) - F{t, u{t)) 9 0 

^ u{t) € C{t) 

, 'w(O) = uo , 

with an initial data uq € C'(O) and F : / x R ^ R a set-valued map taking nonempty com¬ 
pact values. Here, N(C'(t), n(t)) stands for the proximal normal cone to C{t) at the point 
u{t). This differential inclusion can be thought as following: the point u{t), submitted to 
the field F{t,u{t)), has to live in the set C{t) and so follows its time-evolution. 

Let us first give a brief overview of the study for such problems. The sweeping processes 
have been introduced by J.J. Moreau in 70’s (see [24]). He considered the following 
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problem: a point u{t) has to be inside a moving convex set C{t) included in a Hilbert 
space. When this point is catched-up by the boundary of C{t), it moves in the opposite 
of the outward normal direction of the boundary, as if it was pushed by the physical 
boundary in order to stay inside the convex set C{t). Then the position u{t) of this point 
is described by the following differential inclusion 

-u{t) € dlc(t){u{t)). ( 2 ) 

Here we write die for the subdifferential of the characteristic function of a convex set C. 
In this work, the sets C{t) are assumed to be convex and so dlc(t) is a maximal monotone 
operator depending on time. To solve this problem, J.J. Moreau brings a new important 
idea in proposing a catching-up algorithm. To prove the existence of solutions, he builds 
discretized solutions in dividing the time interval I into sub-intervals where the convex set 
C has a little variation. Then by compactness arguments, he shows that a limit mapping 
can be constructed (when the length of subintervals tends to 0) which satisfies the desired 
differential inclusion. 

Thus it was the first result concerning sweeping process (with no perturbation T = 0) 
because dlc(t){x) = N(C(f),x) in case of convex sets C{t). 

Since then, important improvements have been developed by weakening the assumptions 
in order to obtain the most general result of existence for sweeping process. There are 
several directions: one can want to add a perturbation F as written in (1), one may require 
a weaker assumption than the convexity of the sets C{t), one would like to obtain results 
in Banach spaces (and not only in Hilbert spaces). 

In [26], M. Valadier dealt with sweeping process by sets C{t) = \ mt{K(t)) where K{t) 

are closed and convex sets. Then in [10], C. Castaing, T.X. Diic Ha and M. Valadier have 
studied the perturbed problem in finite dimension with convex sets C{t) (or complements 
of convex sets). In this framework, they proved existence of solutions for (1) with a 
convex compact valued perturbation F and a Lipschitzean multifunction C. Later in [11], 
C. Castaing and M.D.P. Monteiro Marques have considered similar problems in assuming 
upper semicontinuity for F and a “linear compact growth”: 

F{t, x) C /3(f)(l + |xl)H(0,l) , V(f, x)€l X M". (3) 

Moreover the set-valued map C was supposed to be Hausdorff continuous and satisfying 
an “interior ball condition”: 

3r>0, B{0,r)cC{t),yteI. (4) 

Then the main concept, which appeared to get around the convexity of sets C{t), is the 
notion of “uniform prox-regularity”. This property is very well-adapted to the resolution 
of (1): a set C is said to be i]-prox-regular if the projection on C is single-valued and 
continuous at any point whose the distance to C is smaller than r/. 

Numerous works have been devoted to study sweeping processes under this assumption of 
prox-regularity. The case without perturbation {F = 0) was hrstly treated by G. Colombo, 
V.V. Goncharov in [12], by H. Benabdellah in [2] and later by L. Thibault in [25] and by 
G. Colombo, M.D.P. Monteiro Marques in [14]. In [25], the considered problem is 

( -du e'N{C{t),u{t)) 

\ u(0) = uo , 

where du is the differential measure of u. The existence and uniqueness of solutions of (5) 
are proved with similar assumptions as previously. 

In an infinite dimensional Hilbert space, the perturbed problem is studied by M. Bounkhel, 
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J.F. Edmond and L. Thibault in [9, 25, 15, 16]. For example in [16], the authors show 
the well-posedness of 

f —du G N(C(t),u(t)}-h F(t,u(t))dt 
\ ( 6 ) 
[ u(0) = uo , 

with a set-valued map C taking r/-prox regular values (for some rj > 0) such that 

\dc{t){y) - dc{s){y)\ < y{]s,t]), yy € h, y s,t € i, s<t (7) 

where /r is a nonnegative measure satisfying 

sup//({s}) < (8) 

The proof uses the algorithm developed by J.J. Moreau with additional arguments to deal 
with the perturbation F and the prox-regularity assumption. 


Indeed the main difficulty of this problem is the weak smoothness of the proximal normal 
cone. For a hxed closed subset C, the set-valued map x —>■ N(C', x) is not in general upper 
semicontinuous and the uniform prox-regularity assumption permits us to obtain this 
required smoothness. To finish we mention that the sweeping processes are also studied 
in the framework of Banach spaces (see [3] and [5]). 

1.2. Sweeping processes on Riemannian manifolds. — As previously explained, 
sweeping processes have been extensively studied in the context of a linear ambient space 
(Euclidean space, Hilbert space, Banach space). However, the proof relies on the catching- 
up algorithm (which only requires a projection operator) and the notion of uniform prox- 
regularity (which is equivalent of an hypomonotonicity property for the proximal normal 
cone). These observations make us want to extend the study of such differential inclusions 
in the framework of Riemannian manifolds, where the metric allows us to define a projec¬ 
tion and where the Hilbertian structure of the tangent spaces is necessary to have a kind 
of hypomonotonicity property. 

In the Hilbertian situation, the finite dimension is not useful and so we aim to work on a 
(possibly infinite dimensional) Riemannian Hilbert manifold (see Definition 2.1 for details), 
which is locally isomorphic to a separable Hilbert space. Let M be such a manifold, we 
consider a sweeping process on M, given by a set-valued map C'(-) and a single-valued 
perturbation / 

N(C'(t), x(t)) - fit, xit)) 9 0 

^ x{t) G C{t) (9) 

^ x(0) = xo , 

with an initial data xq G C'(O). Note that the hrst differential inclusion takes place in the 
fiber-tangent space To study the well-posedness of this problem, we will generalize 

the notion of a prox-regular set (Dehnition 3.7) which will be shown to be equivalent to 
a kind of hypomonotonicity property (see Theorem 3.15). Some difficulties arise from 
the framework of manifolds. Indeed in a Hilbertian space H, a nonempty closed set 
C C H is said r/-prox-regular or prox-regular with constant y if for every point 2 : satisfying 
dc{z) = \nid{z,c) < y there exists a unique point c^ G C such that d{cz,z) = dc{z) and 

cEC 

the map 2 1 —>■ is continuous. However in infinite-dimensional Riemannian manifolds, 
geodesics minimizing the distance between two given points do not generally exist. Thus 
we will make some assumptions so that they do always exist locally. That’s why we will 
give and use a definition of a locally prox-regular set which will take into account the 
geometry of the ambient manifold even if a definition of uniformly prox-regular sets will 
be given in Remark 3.8 in specific manifolds. 
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In addition we will recover an important property of prox-regular sets. Indeed we will 
prove that the metric projection onto a locally prox-regular set is single-valued in its 
neighborhood. 

The notions of proximal normal cones and prox-regularity at a point are already 
introduced in [19] in the setting of a Riemannian manifold. However we would like to 
propose another definition of the proximal normal cone (related to the metric projection 
and better adapted to deal with problems of sweeping processes). Obviously we will show 
their equivalence. Moreover the prox-regularity at a point is not sufficient to treat our 
differential inclusions : we need a ’’larger” concept of prox-regularity which we will call 
local prox-regularity as already specified. Once again these different notions are compati¬ 
ble. Indeed we will prove that a locally prox-regular set in our sense is prox-regular at 
any point in the sense of [19] . 

Finally we prove the main result: 

Theorem 1.1. — Under some assumptions about the Riemannian Hilbert manifold M 
(Assumption (2.9)), we consider a Lipschitz bounded mapping f and a Lipschitz set-valued 
map C{-) taking nonempty and locally prox-regular values (Assumption 4A). Then for all 
xq G 0(0), the system (9) has a unique solution x G M). 

To our knowledge, sweeping processes on (finite- or infinite-dimensional) manifolds have 
not been previously studied. However other differential inclusions on Riemannian mani¬ 
folds have already been treated. For example, in [18], the author is interested in control¬ 
lability properties of autonomous differential inclusions: 

x{t) G F{x{t)) 

in compact finite-dimensional Riemannian manifolds. The map F is supposed to be Lip¬ 
schitz continuous and takes nonempty convex and compact values. Then in [13], the 
authors prove that there exists at least one solution of the problem: 

x{t) G F{t,x{t)), 

with an initial condition in a finite-dimensional Riemannian manifold. In this work, F 
is an integrable bounded field, taking nonempty and closed values. In addition this map 
is supposed to be lower semicontinuous or to satisfy the Caratheodory conditions. Note 
that compactness arguments are used in their proof so the finite-dimensional assumption 
is crucial and by such a way, it is not possible to get uniqueness results. 

1.3. Motivation from unilaterally constrained problems. — Let us briefly de¬ 
scribe one of our motivations, coming from problems involving unilateral constraints. 

Even in the Euclidean space several applications for modelling are concerned with 
sweeping processes involving an admissible set 

C{t) := |x G M'’*, Vi G T, gi{t,x) > o| , 

where {gi)i£X is a finite collection of inequality constraints. 

To apply the existence results to the sweeping process with this moving set, we have 
to check that C{t) is uniformly prox-regular (with a time-independent constant g), which 
may be very difficult. However a criterion was given in [27, Proposition 2.9] and was 
used in several recent works ([6, 7, 4]). It would require to check a kind of Mangasarian- 
Eromowitz constraint qualification (MFCQ) of the active gradients (we recall that the 
active gradients correspond to mappings gj satisfying gj{t,x) = 0). Let us specify that 
this condition describes a positive-linearly independence of the active gradients. Of course, 
if the admissible set is also given by equality constraints, 

C[t) := |x G M'’*, Vi G Xi gi(t,x) = 0 and Vi G I 2 , 9 i{t,x) > o| , 
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the same criterion cannot be verified by turning equality constraints into two inequality 
ones (that is to say by writing for i G Xi, gi(t,x) = 0 as gi{t,x) > 0 and —gi{t,x) > 0). 
Indeed a gradient and its opposite are clearly not positively independent. 

So it is important to have another approach for dealing with such admissible sets (in¬ 
volving inequalities and equalities constraints). In the case where the equality constraints 
are time-independent, the new point of view is to consider the manifold given by the 
equality constraints 

M := {x G Vz G Xi, gi{x) = 0}, 

and the admissible set becomes C{t) := {x G M, Vi G l 2 , 9 i{t,x) > 0}. 

Under some conditions, it is doable to check that the manifold M is a smooth sub¬ 
manifold of and so a finite dimensional Riemannian manifold. Consequently it suffices 
to check that the admissible set C'(-) that is only defined with inequalities, is locally 
prox-regular. 

The paper is organized as follows. Section 2 is first devoted to some definitions and 
reminders about Riemannian geometry. Then in Section 3, we extend some notions of 
convex analysis in the setting of a Riemannian Hilbert manifold (proximal normal cone, 
prox-regular set, ...) and we focus on establishing a characterisation of the prox-regularity 
property in terms of hypomonotonicity for the proximal normal cone (see Theorem 3.15). 
Then we prove Theorem 1.1 in Section 4, using the catching-up algorithm and studying 
the convergence of discretized solutions. 


2. Preliminaries and definitions, related to Hilbert manifolds 

For a complete survey on Riemannian geometry and the extension to infinite dimensional 
framework, we refer the reader to [21, 22]. Let us just recall the main concepts. 

Definition 2.1 (Hilbert manifold). — M is a Hilbert manifold if there exists a sep¬ 
arable Hilbert space H and a smooth atlas, defining local charts from M to H. Then it 
follows that for all x G M, the tangent space T^M has a Hilbert structure. Its inner prod¬ 
uct is denoted by {,), {,)x or {, )t„:M cmd the associated norm is represented by \ \ or 
I Such a manifold is called a Riemannian Hilbert manifold if it is endowed with the 

following Riemannian metric: for two different points x,y ^ M 

d{x,y) := inf [ 7 ( 5 ) jds : j is a -curve joining x = 7 (a) and y = 7 ( 6 ) 

The bundle tangent space is denoted by TM and we recall that 

TM = U {x} X TxM. 

xeM 

We denote by V* the Levi-Civitd connection that is the unique connection consistent with 
the metric d. 

In the sequel, M will be a connected Riemannian Hilbert manifold and V* its natural 
Levi-Civita connection. Using this connection, we can define the parallel transport: 

Definition 2.2 (Parallel transport along a curve). — For a -curve 7 : [o,r] ^ 
M and v G the parallel transport of v along 7 is the unique vector field t e^V{t) 

satisfying 



V{t) G Xy( 7 M and 


V*.^7)U(t) = 0 
U(0) = V. 
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Definition 2.3 (Geodesics). — The geodesics are defined to be (locally) the shortest 
path between two points. As a consequence, a geodesic 7 is also a curve such that parallel 
transport along it preserves the tangent vector to the curve, so 

v;7 = 0 . 

For any point x G M and for any vector v € T^M, there exists a unique geodesic 
lx,v ■ I ^ M such that 

f 7x,v(0) =x 
{ 7x,v(0) = V, 

where / is a maximal open interval in M containing 0 (depending on x and v). 

Definition 2.4- — The Riemannian Hilbert manifold M is said geodesically complete if 
every maximal geodesic is defined on M (i.e. I = M). 

From now on, we assume that the connected Riemannian Hilbert manifold M is geodesi¬ 
cally complete. 

Definition 2.5 (Exponential map). — If M is geodesically complete, we recall that 
the exponential map at x, denoted by exp,j,, is defined in the whole tangent space T^M and 
maps each vector v of the tangent space T^M to the end of the geodesic segment [x, y] in 
M issuing from x in the direction v and having length d{x,y) = |n|. 

Definition 2.6 (Injectivity radius). — We define the injectivity radius of M at x, 
ijviix) as the supremum of the numbers r > 0 of all balls B{0x,r) in T^M for which exp^, 
is a diffeomorphism from B{0x,r) onto B[x,r). Then, we denote the global injectivity 
radius of M by 

i{M) := inf iuix). 

x(^M 

Definition 2.7 (Convexity radius). — The convexity radius cm{x) of a point x € M 
is the supremum of the numbers r > 0 such that the ball B{x,r) is convex (in the sense 
that given two points p,q G B{x,r), there exists a unique geodesic in B{x,r) joining them 
and of length d{p, q)). We denote the global convexity radius by 

c{M) := inf cm{x). 

x&M 

Similarly we dehne for a bounded subset B of M, 

i(B) := inf iiuix) and ciB) := inf cm(x). 

x&B x&B 

Note that for x G B, the open ball B{x,c{B)) is not in general included in B. 

Let B a bounded subset of M and x, y G B satisfying d{x, y) < c{B), we define G T^M 
as the unique vector w gT^M satisfying 

7x,w((f^ — U- 

We can note that its norm is equal to the distance between x and y. 

\I^x,y\= d{x,y). (10) 

According to Definition 2.5, we have y = exp 3 ,(F 3 ;^y). In other words, the vector T^^y could 
be also expressed as follows: 

Vx^y = exp"^(y). 

Definition 2.8 (Sectional curvature). — Let x G M and cJx C TxM a two- 
dimensional plane, the sectional curvature K[(Tx) is the Gaussian curvature of the 
section of M whose tangent space at x is ax (in other words, it is the Gaussian curvature 
of the range of Ox under the exponential map at x). 
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Since the interesting notions have been recalled, we can now specify the hypotheses 
about the manifold. 

Assumption 2.9. — We assume that M is a connected Riemannian Hilbert manifold 
which is metrically complete and geodesically complete and M has a locally bounded ge¬ 
ometry : for every bounded subset B C M , 

— i{B) > 0 and c(B) > 0 

— the sectional curvatures are bounded in B: \K{B)\ < +oo 

K{B) := sup K{(Jx) 

X & B 
cTx C T^M 

Thus we define for a bounded subset B C M, 

p{B) :=mhi{i{B),c{B),-K/{2^J\K{B)\)}. (11) 

Moreover the exponential map exp is supposed to be locally smooth onTM and the inverse 
map Dexp~^ : TM ^ T^M is supposed to be locally Lipschitz: for every bounded subset 
B <Z M and for k G {0,1, 2}, there exist Ce{B) > 0 and L{B) > 0 such that 

l|exp||c;=(B) < (12) 

and for all x,y,z G B satisfying d{x,y) < p{B) and d{x,z) < p{B), for all v G TyM,w G 
T,M 

|exp"Hy) - exp-^(z)|r,M < Ce{B)d{y,z), 

\Dexpf;^{y)[v] - D e:i^pf:^{z)[w]\T,M < Ce{B)dTM{{y,v), {z,w)), 
where dxM represents the distance on the tangent bundle: 

druHy, v), {z, w)) := d{y, z) + \Ly^^v - w\t,m 
and Ly^z is the tangential parallel transport from y to z (see Definition 2.11). 


Remark 2.10. — In the finite dimensional case, it suffices to assume that the connected 
Riemannian manifold M is geodesically complete because in that case Hopf-Rinow’s Theo¬ 
rem asserts that M is automatically metrically complete and the closed and bounded subsets 
of M are compact. Since the functions iM, cm clcc continuous, the quantities i{B) and 
c{B) are necessarily strictly positive. In the same way, \K{B)\ < +oo is also satisfied. 

We aim to finish this section with the following definition: 

Definition 2.11 (Local parallel transport). — Let B C M be a bounded subset and 
let x,y G B two different points satisfying d{x, y) < p{B) < c{B), then there exists a unique 
geodesic curve 7 with 7 ( 0 ) = x and 7 ( 1 ) = y and so T^^y is well-defined. Moreover, we 
denote by Lx^y : TxM —>■ TyM the tangential parallel transport defined as follows: for all 
V G TxM, Lx^y{v) = 1^(1) where V is the unique vector field satisfying 1^(0) = v and 
V*.^ 7 )l/(t) = 0. The map Lx^y also defines a linear isometry between TxM and TyM 
with Lff^y = Ly^x- Moreover it satisfies 


Lx^y (ra;,y) — ^y,x- 


( 13 ) 
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3. Technical preliminaries about convex analysis and the squared distance 

In this section, we aim to extend the notion of prox-regular subset to the setting of Rie- 
mannian Hilbert manifold and to understand how the main property of hypomonotonicity 
(of the proximal normal cone) is modified (see Theorem 3.15). 

In all the sequel we consider a manifold M which satisfies Assumption 2.9. We first 
recall some properties about the smoothness of the squared distance: 


Proposition 3.1 (Proposition 2.2 [1]). — For every bounded subset B of M, there ex¬ 
ists > 0 such that for all p,q & B satisfying d{p, q) < p{B): 

\\Hxd^{p,q)\\p <Cp(^tS), ( 14 ) 


where 

with 


Hxd^{p,q) := D'^cfip) 

(p{p) '■=d^{p,q) and ||^||p := sup |'I'[u,u]|. 

v^TpM^\v\=l 


Lemma 3.2. — For all bounded set B included in M, for all r G [0,1] and x,y ^ B 
satisfying d{x,y) < p{B), we have 

< L{B)\Tx,y\^ = L{B)d{x,yf, (15) 

where 'yx,T^,y{r) = expy.{TTx,y) = exp,^(r exp-i(?/)). 

Proof. — Let x,y £ B and consider for r G [0,1] 

lx,r,,yir) = expy,{TTx,y). 

By differentiating in r, we get 

3x,rj:^y{T) = D ex.py.{TTx^y)\Fx,y,Fx^y\. 

Then the assumption (12) gives the expected result. □ 

As a particular case of [20, Theorem 1.2] (with a constant function /), we recall the 
link between the exponential map and the gradient of the squared distance. 

Lemma 3.3. — For every bounded set B C M, for all x,y £ B satisfying d{x,y) < p{B), 
we have 

Vxd^{x,y) = 2d{x,y)Vxd{x,y) = -2rx,y = -2expf^{y). (16) 

For the sake of readability, we give a short proof. 

Proof. — By definition of p{B), there exists a unique geodesic 'yx,rx,y such that 7 x,rj;,j,(l) = 
y. For all t G [0,1), the constant velocity along the geodesic implies that 

dhx,r^,y{t),y) = d{x,y) - t\Tx,y\ = (1 - t)d{x,y). 

By differentiating at t = 0, we get 

{'^xdi3x,r^,y{0),y),jx,r^,y{0)) = {'^xd{x,y),Tx,y) = -d{x,y). 

From \S/xd{x,y)\ = 1 and [Faj^pl = d{x,y), we deduce (16). □ 

Remark 3.4- — With the notations of Proposition 3.1, Vxd‘^{p,q) = V(l){p). We deduce 
from this proposition and from (16) that the map p i->- exp“^(( 7 ) is smooth on B. 

Remark 3.5. — In the same way, 

Vyd^{x,y) = 2d{x,y)V yd{x,y) = -2Ty^x = -2exp“^(3:). 
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Definition 3.6. — Let C be a closed subset of M and x € C. The proximal normal cone 
to C at X is defined as follows: 




Note that N(C', x) is a cone because for all a > 0, for all v G N(C', x), av belongs to 
N(C',x) due to the equality 'yx,v{£) = 3x,av (f)- 

Definition 3.7. — Consider a nonempty closed subset C C M, C is said locally prox- 
regular if for all bounded set B C M, there exists 0 < r]Q < p{B) such that for every 
t G [0, ryfi), X ^ C r\B and v G N(C', x) \ {0} 



Remark 3.8. — Note that we can define a uniformly prox-regular set : a closed subset 
C C M will be called rj-prox-regular or uniformly prox-regular with constant rj if for all 
t G [0, r/), X G C and v G N(C', x) \ {0}, 



However this definition is meaningful if i{M) > 0, c(M) > 0 and \K{M)\ < +oo (for 
example il the framework of compact manifolds). 

Now we want to prove that the local prox-regularity provides a kind of hypomonotonicity 
property. This result is the purpose of the following lemma. 

Lemma 3.9. — Let C be a closed subset of M and assume that C is locally prox-regular. 
Then for all bounded subset B of M, there exists a constant Eq > 0 such that for all 
x,y € C n B satisfying d{x, y) < p{B) and v G N(C, x) 


{v,Tx,y)T^M < Ei3\v\d{x,y)‘^. 


Proof. — Let x,y G C (1 B and v G N(C', x). We can assume |x| = 1 without loss of 
generality. We define B := {m G M, d{m, B) < pb] (where pB is introduced in Definition 
3.7) and we set r := m.m{pB., p{B)). If d{x,y) > rl2, the result is easily proved by (10) 
with Eb = 2lr. 

Consider also x,y with d{x,y) < r/2. Since C is locally prox-regular, for all t G (0,r/2) 




t^ = d{x,px,v{t)f < d{y,jx,v{t)f- 

By a second order expansion and setting m := Px,v{t), we get thanks to (16) 




d{x, m)^ -|- 2d(x, m){'Vxd{x, m),rx,y) 

r^rid^.o, ... .1 
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Consequently, 

-2d{x,m){Vxd{x,m,T:c^y) < 


1 rsT^2 


10 JO L 
fi rs 


-^d 


drds 


- |s=T 


< 



0 JQ 


[Hxd^ {lx,T^,y (r), m) [^x,T^,y (r), 7 a;,r,.^ (r)] 

+ {'^xd?{-lx,T^^y (r), m),(r))] drds. 

Since d{x,'yx,T^,y{T)) < d{x,y) < rjl and d{x,'^x,vit)) <t< rjl, it is clear that 7 x,r,,.j, (t) 
and 'yx,vit) belong to B{x,rj2) C B. Moreover 

'^( 7 x,r,,j,(r),m) = d( 7 x,r,,^ (r), 7 x,^(t)) < r /2 + r /2 <r< p{B). (17) 

Thus by Proposition 3.1, 

|77xa!^(7cc,r,.,(r),m)[7^,r,.,(T),7x,r,.,(r)]| < C. 

Furthermore Lemma 3.2 with (17) gives 

(V,^,(i2(7a;,r,,j,(r),m),7,,,r,,j,(T)) < 2d{-i^^T^,y{T),m)L{B)\T^^y\^ < 2rL{B)\Vy,^y\^. 

Thus by adding the previous inequalities we obtain 

-2d{x,m){Vxd{x,m),r^^y) < Cp(^^^\Ta;^y\'^ + 2rL{B)\Ty;^yf. 

So we deduce from Lemma 3.3 and (10) that 


7 (B) + 2 >-i(e) 


(exp„'( 7 i,»(<)).Ti, 5 ) < 

As exp“^( 72 ;_.i,(t)) = tv, it comes 

CpiB)+ 


d{x,yf. 


{v, ^x,y) < 


2t 


d{x,yf. 


This inequality holds for all t G (0, r/2) which ends the proof by setting Ejs : = 
max(C^(g)/r + 2L{B), 2/r). 

□ 


Remark 3.10. — Note that the constants Eq depend on the geometrical features of M 
but also on C through the constant 

Remark 3.11. — Let C he a closed subset of M and x G C. If v £ N(C, x) then there 
exists a > 0 and A > 0 such that for all y £ C Ci B(x, a) 

{v,r^^y)T,,M < A\v\d{x,y)^. 

Indeed, assume that |n| = 1 (without loss of generality) since v £ N(C', x), there exists 
e > 0, such that x £ Pc(7x,v(£)) by Definition 3.6. Now it suffices to fix B = B{x,e) and 
take a := min(e, p(;B)/2). Thus for all t < a, x £ Pci7x,v{t))- By following the previous 

proof, we can choose A := + 2L{B). Here a and A depend on v and may degenerate. 

The important property of a locally prox-regular set is that these two quantities may be 
assumed to be independent on v, as proved in Lemma 3.9. 

The following statement will be useful in Lemma 3.14 to check that the hypomonotonic- 
ity property (which is detailed in Lemma 3.9) implies the local prox-regularity. Moreover 
its corollary will allow us to show that the projection onto a locally prox-regular set is 
single-valued in its neighborhood. 

Lemma 3.12 ([20, Theorem 1.2]). — Consider a manifold M satisfying Assumption 2.9 
and B a bounded subset of M. 
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— If K{B) < 0 then for every x £ B and along any geodesic 7 taking values in B{x, r) C 
B with r < c{B), 

^d^{-f{t),z) > 2|7(i)|2. 

— //0 < K{B) < 5, with S{t) := tcot(t) and r £ ]t),p{B)], we have for all x £ B and 
along any geodesic 7 taking values in B{x,r) C B, 

> 2S{r6^/'^)\f{t)f. 

Corollary 3.13. — Let B a bounded subset of M containing x, zi and Z 2 . Suppose that 

— zi, Z 2 £ B{x,r) C B with r < c{B), if K{B) < 0; 

— zi,Z 2 £ B{x,r) C B with r < p{B) if 0 < K{B) < 5 . 
then there exists A := A(B) > 0 such that 

(^22,0: -^' 2 : 1 —>- 2:2 r^2,zi) ^ Ad(^Zi,Z2') ■ 


Proof. — We set (p{z) := d'^{z,x). Then with a second order expansion of the function 
and (16), we have 


ip{zi) - ip{z2) 


{V^iz2),exp^^ (zi)) + 


f 


dt‘^ 




where 7 (t) is the geodesic between 7 ( 0 ) = Z 2 and 7 ( 1 ) = z\. With A = 2 (in the case of 
negative curvature) or ^4 = 25(r(5^/^) (if the curvature is at most <5), Lemma 3.12 implies 

ip{zi) -<piz 2 ) > {\/ip{z 2 ),expf^{zi)) +Ad{zi,Z 2 f, 


which can be written as 


(f{zi) - ip{z 2 ) > -2{rz2,x,Lz2,zi) + Ad{zi,Z 2 f. 

By symmetry, we have 

ip{z2) - ipizi) > -2{T+ Ad{zi,Z2f‘. 

By summing these two inequalities, it comes 

(fjZ2 ,z\)Tz.2,^ z^^Z2)Tz^M h Ad{^Z\,Z2^ ■ 

Using parallel transport along the geodesic from zi to Z 2 , 

(r21,X) r2i,22)T2jM = ~{Lzi^Z 2^ ZljXtL Z2,Zi)Tz2M ■> 

which is the expected result. □ 


Lemma 3.14- — Let C be a closed subset of M . Assume that for every bounded subset B 
of M, there exists a eonstant kb > 0 such that for all x,y £ CCiB satisfying d{x,y) < p{B) 
and every v £ N(C,x) 

(v,rx,y}nM < KBlvld(x,yf. (18) 

Then there exists 6^ g]0,/?(B)[ such that for all t < 6^ and for all v £ N(C', x) \ {0}, 

xGPc (7x,^(i)). 

In other words, C is locally prox-regular. 

Proof. — Let B be a bounded subset of M, let us fix x € CnB and v £ N(C', x). Without 
loss of generality assume that |u| = 1. We define B := {m £ M,d{m,B) < p{B)} and 
we recall that p{B) < p{B). Moreover we set 0^ := min if dI{B) < 0 and 

0^ := min arctan(^^)^ if 0 < K{B) < S. We are going to check that for 

t < < p(B)/3 


dc(7x,v(t)) = t. 


( 19 ) 
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We assume that (19) does not hold and then there exists z G C {z x) satisfying 

d{z,'yx,vit)f < 

and in using a second order expansion, we get 

> d{z,'yx,vit)f = + {Vxcfix,'yx,v{^))^'^x,z) 

"1 r cP 


+ 



0 JO 


ds"^ 


d (7a;,ra; ,('S),7a:,'u(i)) 


drda. 


- |s=T 


Note that is well-dehned since d{x,z) < d[x,^x,v{t)) + d[z,'-^x,v{t)) < 2t < p{B). As 
a consequence, z € B and 

d^ 


ff 


j^d {ix,r^,As),lx,v{t)) 


drda 


- |s = T 

< {-S/xd^{x,-ix,v{t))^'^x,z) = (2exp^^(7a,,^(t)),r^,^) = {2tv,Tx,z), 
thanks to Lemma 3.3. Thus with assumption (18), it comes 

d^ 


ff 


-^d {1 x,t,,As)^1x,vA)) 


- |s=T 


drda < 2tK^d{x, zf. 


Moreover 


d{'lx,r^,As)^'lx,v{f) ^ d{x,-/x,v{t)) +d(x,7a;,r,,,(s)) < + s|r,j;,2 

< t + d{x, z) 

< 3t < p{B) < c{B). 


If Ar(,B) < 0 everywhere, Lemma 3.12 yields for t < 

^Aix,Vx^Af\ = =2d{x,z) . 


-^d^{-fx,r^,As),jx,v{t)) 


|s = T 


Consequently, we have for t < 

d{x,zf < 2tKj^d{x, zf 

which is impossible for t < Thus for t < 9^, dc{'yx,v{t)) = t so x & Pdlx^vit))- If 
0 < K{d < d, Lemma 3.12 yields for t € 


n 

u, 3 


d^ 


p[g2^ fl^J'x,z{A^lx,v{f) 


|S=T 


> 2S{3t5^/fdix,zf, 


for every r € [0, s] and every s € [0,1]. Hence, we get 


d(x,zf < — -d{x,zf 


2Kg tan(3t(^^/^) 


5(3t<5V2) 

which is impossible as soon as t < 9^. 


□ 


Combining Lemmas 3.9 and 3.14, we obtain the following characterization of locally 
prox-regular subsets: 


Theorem 3.15 (Hypomonotonicity property). — Let C be a closed subset of M. 
Then C is locally prox-regular if and only if for all bounded subset B of M, there exists a 
constant Eq > 0 such that for all x,y € C (1 B with d{x, y) < p{B) and v G N(C', x) 

{v,Tx,y)Tz,M < Ei 3 \v\d{x,yf. 

As mentioned in the introduction, the notions of proximal normal cone and of prox- 
regularity (at a point) have been already appeared in [19]. Let us recall them. 
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Definition 3.16 ([19]). — For C a closed subset of M and x £ C the proximal normal 
cone to C at x is defined by 

with U := exp^{V), where V is an open neighborhood ofOx- 
Moreover the limiting normal cone to C at x is expressed as follows: 

Nq{x) := < lim Vi : Vi € NQ{xi) and Xj — >■ x with Xj € C > . 

1^2—)^CX) J 

Definition 3.17. — [19, Definition 3.3] A closed subset C <Z M is said to be prox- 
regular at x & C if there exist e > 0 and a > 0 such that B{x,e) is convex and for every 
X G C n B{x, e), ?/ G C n B{x, e) and v G N^{x) with jxj < e we have 

{v,ex.pf^{y)) < '^d{x,yf. 


First we would like to show that the two notions of proximal normal cones are equivalent. 
This result is based on the following proposition. 

Proposition 3.18. — Let C be a closed subset of M and B a bounded subset of M. 
Consider X G CnB and v G T^M. Then v G N(C', x) if and only if there exist 9 G {0,p{B)) 
and A > 0 such that for all y G C H B with d{x, y) < 9 

{v,Tx,y)T^M < Ad{x,yf. (20) 

Proof. — Let v 0 satisfying (20). For all y G C (iB with d{x, y) < 9, we have 

{v,Tx,y)T^M < K\v\d{x,y)‘^, 

where k := ^. Then it suffices us to check that for some small enough t > 0, 

X G Pc(7x,i;(i)), 

which is equivalent to 

dc{lx,v{t)) = t\v\. 

Let assume that it does not hold: dci'yx,v{t)) < t\v\ and choose 2 ; satisfying d{z,^x,v{t)) < 
t|u|. Following the proof of Lemma 3.14, it comes a contradiction for small enough t 
{t < 9^). In conclusion, v G N(C, x). 

Furthermore, if u G N(C', x) then Remark 3.11 yields that (20) is verihed. □ 

Corollary 3.19. — For C a closed subset of M and x G C, then 

iV<?(x) = N(C,x). 


Proof. — Indeed, [19, Lemma 2.1] states that v G N^{x) if and only if there is a > 0 
such that 

(u,exp-^(y)) < |c^(x,y)^ 

for every y in a neighborhood of x. This property is equivalent to (20). We conclude to 
the equality of these cones thanks to Proposition 3.18. □ 

Now let us specify the link between the different definitions of prox-regularity. 

Proposition 3.20. — Let C be a closed subset of M. If C is locally prox-regular (ac¬ 
cording to Definition 3.7 ) then C is prox-regular at any point x G C (in the sense of 
Definition 3.17 ). 
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Proof. — Let x € C and Bq = B{x, 1), then x belongs to the convex ball B = B{x,e), 
with s = min(l,/?(;Bo)/2). Moreover let x,y € C (1 B, it comes 

dix,y) <2e < p{Bo) < p{B), 

since B C Bq. Thanks to Lemma 3.9, for every v G N{C,x), we have 

(r;,exp" 1 ( 1 /)) < Ei 3 \v\d{x,yf. 

Let w G N^{x) satisfying |ic| < e and w = limtCj, where Wi G NQ{xi) with Xi G C and 
Xj —X. For i large enough, Xi G B and |rt;j| < 2e, hence 

{wi,el^pf.l{y))^^ < 2EBed{xi,yf. 

Furthermore, 

(u;i,exp-/(y)),^,. - (w, exp-^(?/)),r 

= {Lxi-,x{wi),L^^^^{expf:l{y)))x - {w,expf^{y))x 
= (Lxi^xiwi) - w,expf;^{y))x 

+ {Lxi^xiwi),Lx,^xiexpf:l{y)) - expf^iy))^ 

0 . 

Indeed the second term tends to zero since x i—> exp“^(y) is smooth on B thanks to Remark 
3.4. Thus by passing to the limit, the previous inequality becomes: 

{w,expf:^{y))x < 2EBsd{x,yf. 

So we conclude that C is prox-regular at x. □ 

Remark 3.21. — In fact, a closed subset C C M is prox-regular at x G C according to 
[19] if and only if it satisfies a hypomonotonicity property in a convex neighborhood of x 
for every vector o/N(C, •) = . 

Moreover, if C C M is prox-regular at x G C, it is proved in [19, Lemma 3.7] that the 
limiting cone, the Clarke cone (that is the closed convex hull of the limiting cone) and the 
proximal normal cone are equal : N(f{x) = = N^{x) = N(C', x) (where the last 

equality comes from Corollary 3.19). 

Then we aim to end this section by proving the following property (which is also well- 
known in the Euclidean space): close enough to a locally prox-regular subset, every point 
has a unique projection onto it. 

Theorem 3.22. — Let C C M be a locally prox-regular set. Then for every bounded 
set B C M satisfying B C C ^ there exists i{B) > 0 such that for any x G B with 
dc{x) < I{B), the projection set is a singleton: 

Pc(x) = {x*}, 

with X* G C. 

Proof. — Let Bi := {m G M,d{m,B) < 1} and (.{B) := m.m.{p{Bi) /2,t{Bi),1) where 
t{Bi) := (we recall that these constants are defined in Corollary 3.13 and in Lemma 

3.9). Consider x G B with dc{x) < I{B), then for every e > 0 satisfying 

dc{x) + e<l{B)<p{Bi)/2 

we can find y G C such that 

d{x,yf < dc{xf + e^. 

Let us consider the map, defined on M by E := d{x,-)'^. Then by Ekeland’s variational 
principle, we can find z := z{e) G C such that 

- d{y,z) < 1, 

- d(x, z)"^ < d{x, yf < dcix)'^ + e^, 
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— and z is the minimizer over C of the function 'll) := F + ed{-, z). 

Note that d{x,z) < dc{x) + e < £{B) < 1 so d{z,B) < 1 and 2 : G Bi. We denote n(x,e) 
the set of points z G C satisfying the three last conditions. Clearly, for ei > € 2 , we have 
n(a;,e 2 ) C n(x,ei). Now consider e > 0 and zi,Z 2 G n(x,e). Then ^ 1,22 G C H Bi and 
by first-order optimality conditions, we are going to prove that for i = 1,2, there exists 
Vi G Tz^M with \vi\T^,M < (for some constant cjs^ > 0) such that 

Wi :=Tz^,x+Vi GN{C,Zi). ( 21 ) 


Step 1 : Proof of (21). 

Let (j) '.= U (Z M ^ H he a local chart around Zi. By definition, we locally have 

Zi := arg min d{x,p)‘^ + ed{p, Zi) 
pec 



arg min d{(p ^iX),(j) ^{P)f + ed{(f) ^{Zi)) 

Pe4>{cnu) 




argmind((/> ^{P)f + ed{(j) ^(P),(/> \Zi)) + I^^c)iP) 

Pe<t>{U) 


where X = (/>(x), P = 4>{p)., Zi = <i>{zi) and I^{c) represents the indicatrix function of 

^{C). 

Since (/>(C') is closed, the characteristic function I^[c) is lower semicontinuous. So we have 
0 G Xd\4>-\X), r\■)){Z^) + ed^{d{ct>-\-), r\Zi))){Z,) + A, 

for some A G X'^^^^{Zi), where and N^ represent the Clarke subdifferential and the 
Clarke normal cone. We refer the reader to [8, Section 3.1.1] for this result. 

As a consequence, there exists Ui G H with \Ui\H < 1 satisfying 

0 = -2[D(j)-\Z,)]*{rz,,x) + eUi + A, 


according to Remark 3.5 . That implies 

[D^-\z,)nrz,,x) - \Pi e Nj(^)(Z0. 

By composing with {[D(l)~^{Zi)\*)~^ we have 

- \{[Dr\Zi)]*)-\Ui) G ([iir'(^.)]*)-'[Nj(c)(^ 0 ]- 

Moreover, 


w G ([Iirn^.)]*)-'[Nj(c)(^.)] ^ [Dcl>-\Zi)nw) G Nj(^)(Z,) 

^ Vr G T^^^c){Zi),{[Dr\Zi)nw),T) < 0 

since the Clarke tangent cone T^^i^^{Zi) to <j){C) at Zi is the polar cone of N^^(^^(Zj). 
a consequence. 


luG (pr'(Z,)]*)-MN^(C)(^.)] ^Vr Gr^^C)(^ 0 ,(«^,^r'(^*)[r]) <0 


w G 


Dr\Z,){T^^C)iZ^)) 


<c 


As 


by definition of a polar cone. Thus, using the change of charts, we obtain that 

'DrHz,){T^^C)iZ^))y = [TSiz,)Y = X<^izi),. 


where we have used the definition of the Clarke tangent cone in [19] and the fact that 
the Clarke normal cone is exactly its polar cone. Since the set is locally prox-regular, we 
know that w G N(C, 2 ,) thanks to Proposition 3.20. We conclude the proof by choosing 










16 


FREDERIC BERNICOT & JULIETTE VENEL 


Vi = — ^{Zi)]*) ^{Ui). In that case, Wi := + Ui G N(C, Zj) and we have 

liiil ^ (the constant depending on cj) and Bi). 

Step 2 : End of the proof. 

Since |iCi| < d{zi,x) + |i;i| < dc{x) + e + < i{B) + Cjg^e, and d{zi, Z 2 ) < d{zi,x) + 

d{z 2 ,x) < p{Bi) from Lemma 3.9, it comes 

{wi,T+ {w2,'B < ‘^Bs^iiiB) + CB^e)d(2:i,^ 2 )^, 

where > 0 . Using the parallel transport, 

hence we conclude to 

{W2 - L2,^22(w^i)>r22,2i)T,2M < 2Ebi{ 1{B) + CBie)d{zi, Z2f. 

By definition of the tCj, it comes 

Z 2 ,X E V 2 Lzi^Z 2 (1^2]^ ,a; “1“ Ul), Ez 2 ,Zi)Tz^M E: 2Ebi {i{B') + CBi d)d[zi, Z 2 ) 
which yields 

(B Z2,X Lzi^Z 2 {rzi,x)-,'Bz2,zi)T,2^ - + CBie)d{zi, Z2f + 2cBied{zi,Z2). 

As 21,22 ^ B{x^i{15)) with £{15) < p{Bi)^ we deduce from Corollary 3.13 that 
Z2^X Zi,X )) Ez2,Z\)Tz 2M > A{Bi)d{zi,Z2f. 

Hence we observe that 


{A{Bi) — 2 EbiI{B) — 2EBiCBie) d{zi, Z 2 ) < 2cBie 
Ife < ^gi^,thenH(Hi)- 2 EB,^(B)- 2 EBiCBie > A{Bi)l2 > 0 since £(B) < r(Hi) 
As a consequence 

^ (A(B,)-2E,AB)-2E,,c,E '■= “■ 

Let 1 0 and Zn G n(x, e„). It then follows from that for m> n, 

d^ZnjZm') E 


AjB,) 

SEgj 


( 22 ) 


hence (zn) is a Cauchy sequence that converges to a point z* € C (M is metrically 
complete). We obtain z* G Ec{x) since d{x,z*) = dc(x}, consequently Pc(a^) 7^ 0- 
Furthermore, taking e = 0 in (22), we can specify that Pc'(3^) = Cl 


Remark 3.23. — Note that the constant i{B) depends on the geometrical properties of M 
hut also on C through the term Eb (which is determined by the quantities rjB see Remark 
3.10 ). 


4. Sweeping process in a Hilbert manifold 

In this section, we consider a Riemannian Hilbert manifold M satisfying Assumption 
2.9. 

Assumption 4-1. — Let C : [0,T] ^ M be a multivalued map satisfying that for all 
t G [0, T], the set C{t) is nonempty and locally prox-regular. More precisely, for all bounded 
set B C M, there exists 0 < pB < p{B) such that for every t G [0,T], s G [0,r/e), x G C(t) 
and V G N{C{t),x) \ {0}, 

X&Pcit) (7a.,^(s)) • 

Furthermore, C is supposed to be a Lipschitz map: for all t,s & [OjP]? dH{C{t),C{s)) < 
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Assumption 4-2. — We assume that f : {t,x) G [0,T] x M i-A f{t,x) G T^M is bounded 
and satisfies a Lipschitz regularity in the following sense: 

f [0, r] X M ^ TM 
^'1 {t,x) H- {x,f{t,x)) 

is Lipschitz with constant Lf. More precisely, for all bounded set B of M, for all (ti,t 2 ) £ 
[0,T]^ and {xi,X 2 ) G verifying d{xi,X 2 ) < p{B) , 

dTM{g{tl,Xi),g{t2,X2)) := d{xi,X2) + \L^^^a^fif{ti,Xi)) - f{t2,X2)\ 

< Lf{\ti - t 2 \ + d{xi,X 2 )). 


Theorem 4-3. — Under Assumptions 2.9, 4-1 o.nd 4-2, for every xq G M, for all T > 0, 
there exists a unique solution x G M) such that 


x{t) + N{C{t),x{t)) 3 f{t, x{t)) for a.e. t £ [0, T] 
x(0) = Xq. 


(23) 


The proof is a mixture of the classical one in a Hilbert space (see the papers cited in the 
introduction) and of the previous results extending arguments in a Riemannian context. 


Proof. — First step: Construction of discrete solutions. 

We set Bq := M U B{xQ,2T\\f\\oa + KlT)- We can choose T < T such that 2T'||/||oo + 
KlT < min{r]isQ/2,£{Bo)), where i{BQ) is defined in Theorem 3.22 by replacing C with 
(7(0). Note that the quantities ^(-Bo) are the same for all the sets C{t),t G [0,T] 
thanks to Remark 3.23 and Assumption 4.1. 

Then we define B := M n R(xo, 2T||/||oo + KlT) C Bq, hence p{B)l2 > p{Bq)/2. As a 
consequence, it yields 

2r||/||oo + KlT < r,Bj2 < p{Bq)/2 < p{B)/2. 

We fix a time-step h = Tjn which obviously satisfies 


h\\f Woo < p{B)/2 (24) 

So we consider a subdivision of the time-interval J = [0, T] defined by tf = ih for i G 
{0, ..,n}. We build (x”)o<i<n as follows: 


Xq = Xq 


Xi+i ^ Pc(t: 


i + l 




First let us check that this scheme is well-defined. 


(25) 


Lemma 4-4- — Bor all i G {0,... n}, x'f and for all i £ {0,... n — 1} ^xv-,f{t'^,xv-){h) are 
well-defined and belong toB. Furthermore for i G {0, ...n —1} , d{xf,x'f_^^) < p{B)l2 and 
d{x2^i,'yxYj(tY^xf){h)) < p{B)l2. 

Proof. — First Xq = xq £ B and since /i||/||oo < p{B), 7a;jf,/(tjj,xg)(^) exists. Moreover 
d{xo,^x^ji^t^^x^){h)) < hWfWoo hence ^x^j{t^,x^){h) £ B. Assume that for some i < n, x'f 
and 'yxv-j{t‘^,x‘^){h) are well-defined and belong to B. Let us show that and (if i < n—1) 
Ix^ ,xv- ){h) satisfy the same properties. Since d{x'f ,^x^j{t‘^,x‘^){h)) < h||/||oo and 

dcit^^filxf) < dH{C{tf),C{tf^^) < KlL then 

^C{t"^u{'lx’?j{t^,x^){h)) < d{x2,^x’?,f{t'?,x^)ih)) -£ dc{t'?_^^){Xi) 

<h{Wf\\oo + KL)<i{Bo). 

In addition Bq fl (^(t^Yi) ¥" 0- Otherwise, for all z G C{t'^_^_i), d{xo,z) > 2T||/||oo + 
KlT and so dc'(t 7 +i)(a^o) > 2r||/||oo + KlT. Yet dc(t-^^)(xo) < ^^(<7(0), (7(t(Yi)) < 
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KlT. Consequently / 0 according Theorem 3.22, so exists. 

Moreover for 0 < j < i, 

d{x'j , < d(3:”,73.nj(^n^3,p(/i)) + d{x'j_^_i,'yxjj(^tj,x^)ih)) 

< HfWoo + 

<h{2\\f\\^+KL)<p{B)/2. 

Thus d{xo,x2^i) < (i + l)/i(2||/||oo +i^L) < T{2\\f\\oo + Kl), hence belongs to B and 
d{x2+i,'yxf,f{tY,xY){h)) = dc{tY_^^)i7xf,f{tY,xY)ih)) < h(||/||oo + Kl) < p{B)l2. 

Since /i||/||oo < p{B)/2, is well-defined. Finally if i < n — 1, 

dixo,lxY+„f{t^,,xY+,)id)) < d{xo,x^^^) + d{x2^^,jx^^^,f{t^^,xY+^){h)) 

< (i + 2)/i(2||/||oo+iPL), 

< nh{2\\f\\^ + KL)< f {2\\f Woo+ Kl), 

so 'y^n f(.n „n \ £ B. Thc proof is also ended. □ 

Now thanks to the points (x”)o<i<n we define two piecewise maps x” and /” on J in 
the following way: 


Vt G Ji := [ih, {i + l)h[, r{t) := /(C, xD G Tx^M, /"(f) := /(t", x") (26) 


and 


Vt G Ji := [f/i, (i + l)/i[, x"'(t) := 


2;n 

i ’ i+1 


t — ih 
h 


,x"(f) :=x;(. 


(27) 


The function x"' is continuous on J = [0,T]. Note that Tx’?,x’i_^^ is well-dehned since 
d(x"+i,xf) < p{B)/2 (according to Lemma 4.4). 

Moreover we define two other piecewise maps r"' and 0”' on J in the following way: 

Vt G Ji := [ih,{i + l)h[, T^{t) :=t^r"(r) := T and 0^(t) :=t”+i,0^(T) := T. 


Second step: Boundedness of the discretized velocities. 

We claim that the discretized velocities are uniformly bounded. Indeed for n, i £ {0,..., n — 
1} and t £ Ji, the scheme (25) with (27) gives 


|x"(f)| = 


h 


' X. 


t — ih 
h 


So, let 


h 


V := sup/i M(x”,_i,x”) < 2||/||oo + Kl < CO. 


^ I — 

_ ^ I I 

<2||/||oo+i^L. 


Thus we observe that 


Vn and Vt, x^(t) £ B. 


(28) 

(29) 

(30) 


Third step: Differential inclusion for the discrete solutions. 

We are now looking for a discrete differential inclusion satisfied by the function x"'. More 
precisely, we want to check that for almost every i £ { 0 , ...,n — 1 }, we have 


‘-i+l’' 


+ /i/(t((j_i,x(Yi) + R(h^) £ N(C(t^_^i),x^_^i), 


(31) 
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with R(/i^) G and |R(/i^)| = 0{h‘^). By definition of the scheme (25), we know that 

(taking e = 1 in the definition of proximal normal vectors and writing 'yxv-j(tv-,xv-){h) = 
JxyA^) where i; = ) 

A+1 )• (32) 

Note that is well-dehned since d(xf+i,j(in^^n)(/i)) < p{B)/2 (accord¬ 

ing to Lemma 4.4). Thanks to the smoothness of the exponential map and of the geodesics, 
let us check the following equality: 


Lemma 4-5 

r 


,n)(h) =^x^.^,x^ + hDe^pJ {xA[f{A,xf)]+Ri{h^), 


(33) 


Proof. — By setting yf := '^xv-j[tV',x‘^){h)-: we recall that 

Now we define for s G [0, /i], 7 ( 5 ) := exp^n x”)), so by using a first order expansion, 

it comes 


exp J , iVi) = exp J (xf)+ D exp^i {-f{s))[j{s)]ds. 


"i+i 


“i+i 


i+i 


Furthermore 

rh 


[ DexpJ (7(s))[7(s)](is - hDexpJ ix2)[f {tf, x^] 

Jo 

[ DexpA {j{s))[j{s)]ds- [ Dexp"i ( 7 ( 0 ))[ 7 ( 0 )]ds 
/o Jo 

f'h 

ds 


< 


< 


< 


< 


[ D exp^i ( 7 ( 5 )) [ 7 ('S)]c^'S - D exp^i ( 7 ( 0 )) [ 7 ( 0 )] 

fh 

/ Ce{B)dTM{{'l{s),'y{s)), ( 7 ( 0 ), 7 ( 0 )))ds according to Assumption 2.9 

Jo 

Ce{B) j (d(7(0),7('S)) + l-^^7(s)^7(o)(7('S))-7(0)|)ds 

Jo 

Ce{B) [ {s\f{A,xA\+L{B)d{'y{0),'y{s)f)ds using 
Jo 


Lemma 3.2 


<c,{B)mooh^+Lmf\Ah^). 

i-h 

Thus fixing Ri(/i^) := / Dex.p~} ( 7 (s))[ 7 ('S)]Js — hDex.p~f (x”)[/(t”,x”)], we can 

Jo 

write 


which is the desired result. 


exp^n (yf) = exp . (x^ + /iT>exp . ixf)[fiA,x2)] + Ri(/i ), 

i + l 1+1 1+1 


□ 


With Assumption 2.9, the map JJexp J is Lipschitz and so it comes (using the nota- 

^i +1 

tions introduced in this assumption) 

^exp-i ^(x")[/(t",x")] - L'exp-i ^(x"+i)[/(t"+i,x"+i)] 

< Ce{B)dTM {g{tf,x\),g{tf^^,x'f^A) 

<Ce{B)Lf{h + d{x^,x^^A 

< Ce{B)Lf{l + V)h := Cih, 
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for some numerical constant Ci > 0, where we used Assumption (4.2) (on g introduced 
there) and (29). Thus (33) becomes 


+hl)exp^i^(xr+l)[/(c+l, 
= r,r+n-r + /^/(*m,^r+i) + R(/i')- 


Since for every x € M, D exp^^{x) = [T>exp^(0)] ^ = Idx^ 
deduce (31). 


M- 


<+i)]+R(h2) 

From this and (32), we 


Fourth step: Cauchy sequence. 

Let n,m two integers, by (30) for all t G J, x"'(t) and x™'(t) belong to B and so 
d{x^(t), x^(t)) < d{x^{t),xo) + d{xo,x^{t)) < p{B) j2 + p{B)l2 < p{B). We consider 
t ^ {T”(t), 0”(t), r™'(t), ^"^(t)}. From (31) and Lemma 3.9, we get for y = x™'(t) € B, 


^ ^B\^x^(e^{t)),x^{T^(t)) + hf{9^{t),x^{9^{t)))\d{x'^{9^{t)),x'^{t))‘^ 

+ 0{h^)d{x^{9^{t),x^{t))) + Oih^)d{x^{9^{t),x^{t)f. 

Since 

d{x^{9^{t)),x^{t)) < 1 + d{x^{9^{t)),x^{t)f, (34) 

it comes 


{^X^{e^{t)),X'^{T”{t)) + ^/(^"‘(^)) 3:”'(0"'(t))), F3.n(0n(i))^3,m(t))a,n(0n(4)) 

< Esh{V + ll/lloo + 0{h))d{x^{9^{t)),x^it)f + 0(^2). 
As 

d{x'^{9^{t)),x^{t)) < d(x”(r(t)),x"(t)) +(i(x'^(t),x"*(t)) 

< Fh + (i(x'^(t),x™(t)) 

we have 

d(x’"(r(t)),x'"(t))2 < 2d(x’"(t),x”^(t))2 + 0{h^) 

and as a consequence 

(Fa,n(0n(i))^3,n(rn(t)) + /l/(0” (t) , x” (^^ (t) ) ) , Fj-n (gn (j)) (f) ) (Qn (j)) 

< 2EBh{V + ll/lloo + 0 (h))d(x’^(t),x™(t ))2 + 0{h^). 


(35) 


Let us split the left term in Init) + IIn{t) with 

In{t) '■= (ra;'i( 0 'i(t)),x"(T"h))) 

and 

Unit) ■= h(/(0"'(t), x”(0"'(t))), F3,n(0n(j)) 

By symmetry, with changing the role of x^ and x'", we obtain 

Im{t) + limit) < 2EBh{V + ll/lloo + 0(h))d(x”(t), X™(t))2 + 0{h^). 


Thus 

Init) + Imit) Ellnit) + limit) ^ ^ ^ 0(h))d(x"(t), (t))^ + 0(h). (36) 

Let us study all the quantities In,Im,IIn and Ilm- 
For we use the parallel transport from x"'(0"'(t)) to x”'(t) to obtain 
Init) = (ra;'i(0'i(t)),3;"(T"(t))5 

= —/i(x”(t), [Fa-n (Qn (j)) (j) ] ) ^.n (j) , 


( 37 ) 
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Since 






i))] - 

= -hx^{t). 


(38) 


Indeed the curve x”(t) is exactly the geodesic between x"'(T"'(t)) and x”'(0"'(t)) with ve¬ 
locity x"'(r"'(t)))//i. Similarly for m, we have 

Im{t) = x”^ {9”^ (t)) ,X-^ {t)\) X”^ it) ■ (39) 

Concerning I In and Ilm, we split them in two parts Iln = + Bn with 

:= ^(-ha;"(6("(t))->-x"(t) [f i't), {0^ (t)))] j 

x^ {t) ,x”^ (t)) x^ (t) 

and 

-Bn := [/(^”(^)) ^"‘(^"'(^)))] ) 

The hrst term An is estimated as follows: 

l^nl < ^ll/l|oo|r3,n(j)^3,m(j) “ (6»" p))-lx" (t) [rx** (6»" (t)) ,0;"* (t) ] I T^n (j) M 

< /i||/l|oo|exp"J^^^(x™(t)) - Ta:-(0-p))^x-p)[exp"J(g„(^))(x”*(t))]|r,n(t)M 

< ^^ll/lloo|Vxd^(a:”(t),x™(t)) - T3.n(0n(i))^^n(t)[V^(i2^x’"(r(t)),x"*(t))]|T,„(t)M 


according to Lemma 3.3 


<-/i||/|U sup ||iL,d^(a:-(s),x-(t))|Un(,)d(a:-(r(t)),x-(t)) 
^ SG [!,(?"(*)] 


1 


1 


< -Cp^B)h\\f\\ood{x^i9^{t)),x^{t)) < -C,^B)h^\\f\UV, 
where we used Proposition 3.1 and (29). Similarly for Am, we obtain 

|A^| < lc,^s)h^\\f\\ooV. 

Finally with (36),(37) and (39), we obtain 

{x (^))-hx"(6'*(t))—(i) [rx"'(0"'(t)),x™'(t)] )x'*(t) 

— (i:™'(t), -L3;m(gmp))_^3,mp) [l'{9^ {t)) ,X^ (t)]) X”^ (t) 

\Bn + Bf, 


(40) 


(41) 


< 


h 


+ iEeiV + ll/lloo + Oih))dix^{t),x^{t)y + 0{h) (42) 


It remains us to estimate \Bn + Bml- With a parallel transport, we get 
Bn + B 


(-^x"(6"(l))—lx"(l) [/(^ y^iX {0 (t)))] x'^{t),x'^{t))x"^(t) 

+ {Lx^{9^{t))^x^(t) ,I'x^{t),x^(t))x^{t) 

(-^x’*(t)—lx"*(O-^x"(0"(O )—[f (d^^X (0 (^)))] ) ^ x'^ (t) jX’^ (t)} x'^ (t) 

+ (-^x"*(6»"*(i))-13;’"P) [/(^"*(^)5 ^"*(^™(^)))] 5 

since Lx’^(t)->-x’^(t)i^x’^(t),x’'^(t)) = ^-np). Consequently, 

--= (-^'x"(t)->-x”i(oBx"(e"(0)-i3;"p) [f (t)j x^(0^(t)})] 

~ -^X’"(0’"(t))^X’"p) [/(^™’(^)) ^™’(^™^(^)))] ) 3.np))3,m(i) 

< |-^x'>(t)->-x’"p)-^x"(6»"(0)-la;"(0 lf(^"^y)jX^(0"'(t)})] 

~ -^x’"(0"*(t))->'x’"(i) [f(^^(i}jX^(0"^(t)))] |r2,m(t)M d(x^(t), x^(t)}. 
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We estimate the last norm by making appear intermediate points. Indeed, using the 
Lipschitz regularity on / (Assumption 4.2) and (29), we have 

<Lf{h + d{x^{t),x^{e^{t)))) 

<Lf{l + V)h<h. 


Similarly, 

and 

\Lxr^{t)^x^{t) - f{t,x'^{t))\ < d{x'^{t),x'^{t)). 

So we can conclude that 

\Bn + Bm\ ^ (h +d{x^{t),x^{t)))d{x'^{t),x'^{t)) 

< hd{x'^(t), x"^{t)) + d{x'^{t), x'^{t)Y 

< h{l + d(x™(t), x”(t))^) + d(a;”*(t), x"'(t))^ with(34). 


Finally, (42) becomes 

(X 1 L ^Qn \p x'^ (6'^ it)) ,X'^ {t)\') X'^ (t) 

— (x”*(t), [hj-m (j)) 

< (C '2 + 0(h))d(x”(t),x”^(t))2 + 0(h), (43) 

for some numerical constant C 2 (not depending on n and m). For the first terms, since 
the discretized velocities are uniformly bounded, it follows that 

= exp"i(g,,(^^^(x”'(t)) - exp“i(^^(x”'(t))| 

— p{B)dix'^{9^(t)), x'^{t)) < h (same arguments as in the proof of (40)). 

Hence 


{X (6" (t)) jX™-(t) ] ) X" (t) 

{x if') X"^ {t) ,X'^ {t)^ X"^ {t) ~\~C){}l). 

Producing similar reasoning for the symmetrical quantity, it comes 


(x {t), r {t) ,x'^ {t)) x^^ [t) {x (^)) rx"*(t),x"'(t) )x"*(t) 

< {C 2 + 0{h))d{x^{t),x^{t)f + 0{h), 

which gives (according to Lemma 3.3) 

{x^{t),S7xd'^{x'^{t),x^{t)))x:n(^t) + (i™(i), Vxd2(x”"(t),x’^(t))),^m(4) 

< 2(^2 + 0{h))d{x^{t),x^{t)f + 0{h). 


Since Vxd?{x"^{t), x'^{t)) = Vy(i^(x"'(t), x™(t)), we conclude that 

j^d{x^{t),x^{t)f < 2 ( 6-2 + 0{h))d{x^{t),x^{t)f + 0{h). 

As usual, Gronwall’s Lemma implies that 

sup d{x'^{t),x'^{t)f = 0{h). 
tG[0,T] 


(44) 
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and so is a Cauchy sequence of C'^([0, T], M). Since M is supposed to be metrically 

complete, the sequence (x”)^ also strongly converges to a function x € T], M). 

Furthermore 

sup d(x”'(t),x(t)) = 0(/i2) (45) 

46[0,f] 

Fifth step: The limit function x is solution of (23). 


Let to G [0)^]) there exists ^ : Uq ^ H a chart where ?7o is an bounded open set of M 
containing x(to)- As Uq is an open set, 

3r < p{Uo), U := B{x{to),r) C Uq. 

Obviously, r < p{Uo) < p{U). Since x is continuous on [0,T], there is an open interval 
J C [0,T] containing to such that 

Vt G J, d{x{t),x{to)) < (46) 

Moreover as x”’ uniformly converges to x on [0,T], there exists ho such that for all h < Hq, 
for all t G J 


d(x"'(t),x(t)) <-, Vh <-, Kih < - and K^h < (.{U) (47) 

where (.{U) is defined in Theorem 3.22 by replacing C with C'(to). Note that for all h < ho 
and for all t & J, x”'(t) G U. 

Recall that x” satisfies equation (31) on J, which is 

^ 3- f {0"“(t), x^{0'^{t))) + R(/i) G N(C'(0"'(t)), x"'(6i”(t))). 

By Lemma 3.9 and thanks to the boundedness of / and of the discretized velocities 
(28), there exist constants 8 and En such that for all c” G n U satisfying 

d(c^,x^(0^(t))) < p{U), 

ih ^ X'^ (6'^ [t)) ,X'^ {t" (t)) 3~ f {0 (t))X {0 (t))) T R(h) , r,j,n 

</3Ff;d(x"(r(t)),c")2. 

So using a parallel transport with (38), 

{—xE{t) + [/(^"'(^)j^”(^"'(^)))] +R2(/i), 

< f3Eud{x'^{0'^{t)),c^)^, (48) 

where |R 2 (/i)| = 0{h). Now in order to apply usual arguments in the Hilbertian context, 
we use the chart to work in H. So we define y” = <l>(x"'). By this way, and y” are 
bounded sequences and y” strongly converges to y := <h(x). Then we know that up to a 
subsequence, we can assume that y” *—weakly converges to y in L°°{J,H) which implies 
that y"" weakly converges to y in since J is bounded. Consequently by Mazur’s 

Lemma, there exists a subsequence z” G L^{J,H) satisfying for almost every t £ J 

z'^{t) G Conv (y^{t), k > n') -^ y{t) (49) 

\ / n—>-oo 

where Conv stands for the convex combinations and 

z”- >y in L\J,H). (50) 

n—>-oo 

We now look for proving for a.e. t (when (49) holds) and all c £ C{t) nU satisfying 
d{x{t),c) < r/4: 


(-x(t) + f{t, x(t)), < pEud{x{t),cf. 
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Let h < ho, t £ J (such that (49) holds), c G C{t) fl U satisfying d{x{t),c) < r/4. Thus 
< dHiC{t),C{9^{t))) < K^h < liU) by Assumption 4.1 and (47). Indeed, 
we recall that the quantity i{U) is the same for all the sets C{t),t G [0, T] thanks to 
Remark 3.23 and Assnmption 4.1. Consequently Theorem 3.22 implies that there is G 
Pc( 0 -(t))(c)- Since 

d(c, c^) = dc(B^(t)) (c) < Kih < r/4, (51) 

we deduce from (46) that 

d{c^, x{to)) < d{c^,c) + d{c, x{t)) + d{x{t),x{to)) < 3r/4 < r. 

In other words, c"' belongs to U. Moreover with (47) we obtain 

d(c”, x”(r (t))) < d(c^, c) + d{c, x{t)) + d{x{t),x'^{t)) + d{x^{t),x^{9^{t)) 

< r/2 + r/4: + Vh < r < p{U) (52) 

and 

d{c, x"'(0"(t))) < d{c, x{t)) + d{x{t), x"^{t)) + d(x"(t), x"'(0"(t)) 

< 3r/4 < p{U). (53) 

Firstly, the inequality (48) is satisfied thanks to (52). Secondly, 

|(L2,n(gn(j))^3,n(^) [f{e^{t),x^{e^{t)))]+R2{h) , 

< (||/||cxD + Cl(/l))|La,n(0n(i))_^3,n(i)[r3,n(gn(4)yc'» “ (Qn J | 

< (ll/lloo + 0(h))| exp"/(g,(^^)(c’") - exp"/(g„(^^)(c)| 

< (ll/lloo + 0{h))Ce{U)d{d^ ,c) according to (52), (53) and Assumption 2.9 

< (ll/lloo + 0{h))Ce{U)dH{C{t),C{6'^{t))) = 0{h) thanks to Assumption 4.1. 

Thirdly, using (51), we have 

d{x'^{e^{t)),c) <d{x^{e'^{t)),c^) + d{c,d^) <d{x^{e^{t)),c^) + KLh. 

Hence with (52), 

d(x"(r(t)), cf < d{x^{9^{t)), c^f + 2KLhp{U) + [KLhf. 

The previons ineqnality always holds if c” and c are switched (nsing (53)) and as a conse¬ 
quence 

d(x’"(r (t)), c’")2 - d{x'^{9^{t)), cf = 0{h). 

Hence (48) becomes 

{ — x^{t) + [/(0” (t) , x” (0” (t) ) )] + R 2 (h) , La;'* (t))->-a;" (t) rx** (6»" (t)) ,c) T^-n (j) M 

< PEudix^{9^{t)),cf + 0{h). 

For n —>■ 00, since / is Lipschitz continuous on TM and with (45), we have as previously 
{Lx’^{e^{t))^X"it) [f f) , X"'{9^ (t)))] + R2(/i), —La,n(en(ip_^3,n(f)Fa;n(0n(i)) c)T3,n(t)M 

+ {fit,xit)),Tx(t),c)T^(t)M = 0{h^). 

Consequently, we have for large enough k 

{ X (t), x^{9^(t)),c)T^k^^^M 

< ldEud{x\9\t)),cf - {f{t,x{t)),Tx(^t),c)T^,,,M + 0{h-2), 
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which means 


D^{x (t)) [y (0]) 




< l3Eud{x{t), cf - {f{t, x(t)), ra.(t)_c)r,(t)M + 0(/i2) 
because d{x^{0^{t)), c) < d{x{t),c) +0{h^) and so 


kfa-W-l 


-fit), D<^ix'^it)) 


(8^ (t) ^x^ (8^ {t))x 


H 


< PEud{xit),c)^ - {fit, a:(i)), r^(t)_c)T,(t)M + Oih^). 

In addition 

D^ix’^it)) ^ (^-C'a;'=(0'=(t))-i>x'=h)ra;i:(6l'=(t)),c) = [-^^(^(^)) i^x{t),c) + ■ 

Indeed setting Wf^ d^x^(6^(t))—>x^{tfx^(8^{t)),c^ ^x^{t),c have 

iwk) - [L>^>(x(t))"^]* (w) 

1 :i£ 

= D^ix^it))-^ {Wk -Wk + Wk- L^u-.^^ku)w) 


+ 


D^x^it)) - [D^>(x(t)) ^]*(u;) 


= Oihf, 


by smoothness of the chart <I> and because Wk — Wk = C>(/i), Wk — Lx{t)^x^{t)'^ — Oih^). 
From this equality we deduce that 

{-fit), [D<^ixit))-^]* P Eu dixit), c)'^ - (/(t,x(t)),r^(t)_c)T,(,)M + 0(/i^). 

So applying (49) and then (50) give for almost every time t € J and all c G C/ fl C'(t) 
satisfying (i(x(t),c) < r/4, 

i-yit), [D<l>ixit))-^]*r^(^t)fH < l3Eudixit),c)‘^ - {fit,xit)),T,,(^t),c)T,,^t^M, 
which is equivalent to 

i-xit) + /(t,a;(t)),r^(t)^c)r,(t)M < ldEudixit),c)'^. 

As a consequence, Proposition 3.18 yields that for such time t £ J, 

-xit) + fit, xit)) £ N(C(t), x(t)). (54) 


We have shown that for every to £ J, there is an open interval J where inclusion (54) is 
satisfied almost everywhere. By a compactness argument, J can be covered by a finite 
number of these intervals so finally (54) holds for o.e. t £ J. Furthermore by dividing 
I = [0, T] into small intervals of length | J| = T, we can again follow these five steps and 
obtain the same result for a.e. t £ I. In other words, the function x (so built in the whole 
interval I) is a solution of (23). 


Sixth step: Uniqueness 

Let xi, X 2 be two solutions of (23). Since xi and X 2 belong to W^'°°iI,M), for all 
t £ I, xiit),X 2 it) £ B, where B is a bounded set of M. Furthermore setting 


Pirn 

2max(||xi||oo, ||® 2 ||oo) 


> 0 


we have for i = 1,2 and for all t < a, dixiit),xo) < /?(B)/2 and so dfiit), X 2 it)) < piB). 
It comes from Lemma 3.9 that for a.e. t < a. 


(-xi(t) + /(La:i(t)),F,^j(4)^3,2(i))r,^(,)M < E^lxiit) - /(t,a:i(t))|d(a:i(t),X2(^))^ (55) 
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and 

i-Mt) + < EB\x 2 {t) - f{t,X 2 it))\d{xi{t),X 2 it))'^. (56) 

Moreover with (13), xi(t),x 2 (t)) = -'^X 2 {t),xi{t)- Thus (55) becomes 

{Exi{t)—>-X2{t) (®1 (0 ~ fi^iXl {t))), ^X2{t),xi{t))T^^(^t)^ 

< Eslxiit) - f{t,xi{t))\d{xi{t),X 2 {t))^. (57) 

By summing (56) and (57), we have 

(Txi(t)—>-x2(I) (Tl (0 f {i} Xiit))) (af2(t) f (diX2{t))):^x2{t),xi{t))Tj.^(^^^M 

<EBEd{xi{t),X2{t)f, (58) 

where F = 2||/||oo + ||®i||cxd + 11^2||cx)- Assumption 4.2 implies that 

\f{t,X 2 {t)) - Lx^(t)^x 2 {t){f{i^xi{t)))\ < Lfd{xi{t),X 2 {t)). 

As \rx 2 {t),xi{t) \ = d{xi{t),X 2 it)), we have 

{f{t,X2{t)) - Lxj^(^t)^x2{t){f{t,xi{t))),Tx^{t),xi(t)) < Lfd{xi{t),X2{t)f. (59) 

It follows from (58) and (59) that 

(Txi(t)— >-x2(I) (^'l (^)) ~ ®2(0) Ta,2(t),xi(t))Tj.2(t)M ^ {EbE L f) d{xi(t) , X2{t)) , (60) 

Furthermore 

mid{xi{t),X2{t)f = {xi{t),Vxd?{xi{t),X2{t)))T^^f^^^M 

+ (X2 (t) , Vy (t) ^ a;2 (t) ) (t)M 

— {xi{t), —2T x^{t),X2{t))T^^(t)M 

+ (x2(t), X2{t),xi(t))T^^(-i)M 

‘^{^xi(t)^X2{t){x\{t)) ~ X2{t)X2{t),xi{t))T^^(^)M ■ 

We deduce that 

\^{d{xi{t),X2{t)f < {EbF + Lf)d{xi{t),X2{t)f. 

which is equivalent to 

^{d{xi{t),X 2 {t)f < 2 {EbF + Lf) d{xi{t), X 2 {t)f. 

With the help of Gronwall Lemma, we conclude that xi = X 2 in [0, a] since xi(0) = xq = 
2 : 2 ( 0 ) and obviously this equality holds in [0,T] with the same arguments. □ 
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